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ABSTRACT 


Joint probability distributions and conditional means 
functions are used to analyze turbulent wind and temperature 
data ohtained over natural ocean waves. The data were 
peubected aboard R/V FLEP during BOMEX 1969. 

Stable conditions prevailed during all time periods 
eonsidered. The emphasis is to use the statistical pro- 
cedures to describe the wave's influence on the temperature 
fluctuations at the 8 meter level. 

Results are compared with predictions given by potential 
flow theory. The results consistently demonstrate that waves 
exert an influence on the temperature fluctuations. In 
addition, it is observed that under near neutral thermal 
Peat fication me poomonred periods of stable thermal 
stratification the temperature fluctuations approached those 


expected from potential flow. 
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I. INTRODUCTION 


Turbulent and wave induced temperature fluctuations over 
ocean Waves may be described by theoretical models or inferred 
from observations. These descriptions are necessary to 
develop an understanding of a wide class of applied problems 
which require a knowledge of sensible heat transfer and the 
moture Of temperature fluctuations. 

Vertical turbulent transfers of both sensible and latent 
heat in the boundary layer are influenced by modifications 
induced in the overlying air mass by the surfaces below. 

Most present theories and procedures for estimating these 
transfers are based to a large extent on empirical formulae 
developed from data obtained over Landa livese ane lude 

assumed profiles for the mean wind, temperature and humidity. 

However, recent theoretical and observational studies 
that will be described in the background section suggest that 
transfer processes over ocean waves may be influenced, and 
not necessarily in a simple manner, by the waves themselves. 
Therefore, one cannot assume that phenomena which are fairly 
well understood over land are also understood over water 
until the wave's influences are adequately determined. 

Saewap peed problem in whieh a better description of the 
PoMipemaauuremriuctuattons Over waves is required is the 
Select hacer) propagation problem. These studies require 
a knowledge of the eae idexim order to predict the 


Pearce =Seimtdalation of optical wave propagation. 





Temperature fluctuation would influence this index and 
thereby optical wave propagation. Optical wave propagation 
is presently based on stability parameters which ignore the 
influence of waves. 

The purpose of this study is to examine turbulent features 
‘of the temperature structure over natural ocean waves. The 
observational data were obtained aboard R/V FLIP (Floating 
Instrument Platform) during BOMEX (Barbados Oceanographic and 
Meteorological Experiment) during the period 19 to 28 May 
1969. The data examined are from simultaneous measurements 
of velocity components (u,v,w), temperature (T) at fixed 
levels above mean sea level and wave height (n). 

The analyses methods used in this Study are the same as 
those in an investigation of velocity fluctuations over waves 
by Davidson and Frank (1973). However, temperature fluctua- 
tions are the parameters of interest in this study whereas 
the velocity fluctuations were considered in the referenced 
study. Statistical procedures applied in these analyses are 
joint probability density function and conditional mean 
meer ion. The former will be referred to as (JDF) and the 
Pameerteasme (CME) in the remainder of this thesis. This study 
was the first application of the JDF/CMF analysis to tempera- 
ture data over ocean waves. Results will be discussed in 
relation to observational studies by Volkov (1969), Phelps 
and Pond (1971) and Pond et al. (1971). 

Pmeacteronmmtomtrne JDF/CME statistical procedures this 


ead etnchiuce spectral results. Thompson (1972, S 
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conducted a study of momentum transfer features using the 


same statistical procedures and the same time periods. 
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LL. BACKGROUND 


A study of boundary layer wave induced fluctuations re- 
quires some sort of physical theory connecting observed 
PaOperenesetor those we wish to infer. This physical theory 
may be obtained by developing a model which gives realistic 
results or by observing real data and describing what has 
occurred. The results of Such studies will frequently allow 
inferences of properties other than those that are measurable 
directly by instruments. This background discussion will 
first consider some papers presenting theoretical approaches 


and then review recent attempts to use real data. 


A. THEORETICAL APPROACH 

With respect to theoretical predictions Stewart (1969) 
has Bee cated Bide witha Stable thermal stratification the 
Study of turbulent fluctuations is complicated by the fact 
Midtemotton due Eo both turbulence and internal gravity 
waves is present and the aspects we would like to infer are 
Boeaeivedittferent for each kind of motion. In addition to 
the vexistence of two different types of fluctuations, there 
is a nonlinear coupling between waves and turbulence. 
peewanie (L969), while looking for a tracer of turbulent 
motion, stated that; 


",.. the waves are capable of transporting momentum 


but they are not capable of transporting passive 
scalars per se." 


Therefore, he suggests using a scalar tracer transport <n 


Grider towsee turbulent effects. This in effect implic. that 


M2 





a scalar such as temperature could be predicted by potential 
Fiow theory. Velocity, temperature and heat flux streanm- 
lines shown in Figure 1 are those that would occur under 
Bibapte temperature stratification and potential flow. This 
concept of the coupling between a scalar temperature field 
Over a natural wave field is in contrast to that presented 
by Volkov (1969) which will be reviewed later in this text. 
iimaetneoretteal study, Witting (1972) considered an 
irrotational, plane, periodic, linear progressive wave model 
to study the variation of surface temperature induced by 
Waves and thereby make it possible to estimate the amount of 
phase mart away from potential flow. He concluded that (1) 
when the wave period is sufficiently small(< 3 sec) the 
amplitude of the wave induced temperature fluctuation is 
proportional to the wave amplitude and to the square root of 
the wave period, (2) surface temperature maxima lead wave 
crests by one-eighth of a period when ‘the average surface 
temperature is less than the bulk water temperature and (3) 
surface temperature minima lead wave crests by one-eighth of 
a period when the average surface temperature exceeds that 
of the water below. This theoretical study then recognizes 
Towemunaweeadestrtctnations in the scalar temperature field 


BiGmaaevidat ton Lrom potential flow. 


B. OBSERVATIONAL APPROACH 
Since theoretical studies do not agree, attention is 


ireeoeedmmeroOvanrdrobservational studies. Phelps and Pond (1971) 
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and Pond et al. (1971). examined data obtained aboard R/V 

FLIP during BOMEX 1969. -. However, in their study wave motion 
was not measured and no recognition is given to wave induced 
temperature fluctuations. According to similarity theory 
predictions, the intensities of velocity and temperature 
Meet mations should depend only on the wind speed, the rough- 
ness of the underlying surface and the thermal stratifica- 
tion. In the case of temperature fluctuations, this depend- 
ence can be expressed as 


oT 


“Tt, ee yy) 
where 
* 

tr = (r-T)t/? 
ee men v2 
L=- Tu, /Kkg wT! 
eo von Karman s Constant 
w= w- 7 
w = observed vertical motion 
w = mean vertical motion 
T' = T-T 
IT = observed temperature 
iw = mean temperature 
Z ‘= height above mean sea level 


Therefore, the ratio OT/ZI, should depend only on Z/L. Such 


a dependence has been verified in numerous overland 


LS 





experiments. In particular, the ratio has been observed to 
Tio ede Oormapproximately | in mear néutral conditions. 
Phelps and Pond (1971) observed agreement between this 
oy Pyeresults and the similarity theory predictions with 
Meepect to 2Z/L. Howeyer, they also observed that oT/T, 


Varied as c/u where CG corresponds to the phase speed of 


ke? 
the predominant surface wave. Kitaygorodskiy (1969) and 
Volkov (1969) initially suggested that such a relationship 
should exist. Because oT/T, is inversely proportional to 
the efficiency of the turbulent flow to transport sensible 
heat, the latter findings imply that over a wave surface 
there are considerable temperature fluctuations which are 
not associated with heat transfer. The latter, of course, 
was suggested by Stewart (1969) with respect to wave motion 
Toamieneaesthermally stratified layer. 

The observed relationship between OT/T, and c/u, appears 
in Figure 2. The curves mean square deviations are those 


Gonstructed by Volkov (1969). This suggested curve will be 


compared later to results obtained in this study. 
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20 66 8O 


40 
Cc 
jf Uy (From Volkov 1969) 


Figure 2. on/  . vs us According to Monin Obukhov 
Similarity Theory. 


Phelps and Pond (1971) used data similar to Pond et al. 
(1971) to describe the similarities and differences between 
Eemperature, humidity, and velocity fluctuations. They found 
that over a wide frequency range the temperature waveform 
leads the vertical velocity and humidity waveferm by 20-40 
@eimees, FUE they attributed this to instrument effects and 
Woetbihe Semotion. Again it should be noted that no wave 
firapiewerenavatlable for their eaxperiment. They also con- 
Teded me iniaeenumidity and vertical yelocity are basically in 
Bhaae . No humidity data were available for the present study. 

In another observational study Volkov (1969) reduced 
Nectenumieditertanean data according to similarity theory 


Mio edcrmuoncompare characteristics of turbulence in 
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atmospheric boundary layer above an agitated sea surface with 
measurements in the boundary layer above dry land. He states 
that the main difficulty in interpreting atmospheric measure- 
ments above the sea is the presence of a mobile fluid air- 
water interface which has an irregular variation as a result 
of wind action and airflow-wave feedback processes. The feed- 
back is due to the surface waves disturbing the airflow 
around them and forming waves in the air which interact with 
the velocity gradient of the basic flow. The wave induced 
motion in the air flow causes some energy to be transferred 
from the basic flow to the waves. These processes also 
effect the fluctuation of a scalar such as temperature or 
moisture. In particular Volkov (1969) suggests that an 
observed temperature (T) above a wave surface should be made 
up of- the mean temperature Cie yon EuUMbULen Est luctuations 
Mims y,2,t)) and a imducederluctuidt tone (4 (x,Z)). 
Therefore: 

T = T(z) Sexe i Tt (x.y ,2. t ) 

Huemeromten ot Special Significance in studying the 
structure of turbulence in the atmospheric boundary layer 
above the agitated sea surface are measurements of the fluc~ 
B@ettons wor velocity (u,w), temperature (1) and emphatically, 
wave height (yn). It should be noted that Phelps and Pond 
(1971) and Pond et al. (1971) did not have simultaneous 
measurements of wave ween included in their data. 

By the use of eee analysis Volkov (1969) was able 


momconeimdemtnat the spectra of fluctuations above th¢ 
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agitated sea surface are due to two types. One type of 
fluctuation corresponds to those smaller than the charac- 
Betictwe acales of sea swell and is close to the correspond- 
ing interval of the turbulence spectrum for the air layer 
above land. The other type corresponds to fluctuations of 
the order of the characteristic size of the surface waves 
and is typlifiedby the superposition of wave disturbances on 
ateerandom turbulent field of the fluctuations. The latter 
type was shown in Figure 1 as potential flow. 

Davidson and Safley (1973) in another observational 
Study described the occurrence of micro-thermals, large 
positive temperature fluctuations, and provided arguments 
Supporting the concept of wave influence on micro-thermals. 
A review of the literature associated with wind-wave coupling 
was provided by Davidson and Safley (1972) who strongly 
Suggest that the wave's influence could contribute to micro- 
thermal formation in the case of two proposed mechanisms. 
These mechanisms were the existence of wave-related motion 
in airflow which is too organized to be viewed as turbulent 
MiGdmeinstantaneous distortions of the wind shear over difiier- 
ent phases of the surface wave. Velocity spectral results 
from a Lake Michigan study were used to show that the waves 
foumdumave Dlayed a role in the observed micro-thermal 
phenomena due to wave related vortices in the mean shear. 
However, large thermal instability associated with the 
periods used in the referenced study did not allow a detail 


description of the precise role of the waves. 


ig 





On the basis of 
(1969) and Davidson 
is reason to expect 
influence on all of 


bulence in the near 


BMemresulis of Witting (1972), Volkov 

and Safley (1972) just described, there 
that the wave surface has a considerable 
the statistical characteristics of tur- 


surface layer. Volkov (1969) states this 


fact in the following manner: 


"We therefore succeed in finding some typical features 
and regularities of this influence by the complex 
method of investigating the processes of interaction 
between the atmospheric boundary layer and the agitated 


sea surtrace,. " 


The present study is an attempt to describe the above 


processes through the use of actual data and statistical 


procedures to test concepts put forth by Volkov (1969) whose 


results were. obtained by semi-empirical arguments. 
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TIL. DESCRIPTION AND MEASUREMENT OF DATA 


Pree LOCATION OF R/V FLIP 

The data were from measurements taken by University of 
Michigan personnel during the BOMEX experiment of 1969. The 
data consisted of velocity, temperature, vertical motion 
and wave height measurements at three and eight meter levels. 
The data examined were from 19, 26, and 28 May 1969. During 
this period R/V FLIP was located northeast of Barbados, 
Beegure 3A. 

FLIP's position into the wind was held by a tug attached 
with a bridle and approximately 800 meters of line. The tug 


then steamed slowly downwind. 


B. SENSOR MOUNTING AND INSTRUMENTATION 

ipstrumente Location relative to R/V FLIP is shown in 
Figure 3B. As indicated, all sensors were mounted on the 
vertical mast. 

Temperature fluctuations were measured by a resistance 
Eiemmometer, Figure 4A. The sensor consisted of a 30 ohm, 
0.00015 inch diameter tungsten wire, 0.16 inches long as 
part of a constant-current circuit with a small filament 
current (about 2ma) to serve as a resistance thermometer. 
The frequency response of the system is determined primarily 


by the thermal inertia. 
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Figure 3. Arrangements during BOMEX (A) Geo- 


Baapire position, (B) Sensor location on board 
R/V FLIP. 
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peo ckmararram of the bridge circuit is shown in Figure 


43.7 


An unbalance of the bridge is caused by temperature 
and thereby a sensor resistance change. 

oystem mnolse and sensitivity restricted spectral descrip- 
tion of temperature to frequencies below 20 Hz. Since 
sensible heat flux is associated with frequencies below 10 Hz 
(Deacon 1959), the response was sufficient for this aspect 
of the interpretations. Furthermore, of particular interest 
is the frequency response near the frequency of the dominate 
Surrace wave, approximately 0.1 Hz. 

Velocity measurements were made with a constant-tempera- 
Eure hot-film Macnometer- System using three sensor probes. 
A discussion of this equipment is provided by Portman et al. 
97 G)) . 

Visual estimates of wave height and 10 meter wind speed 
are shown in Figure 5. All periods considered are charac- 


Bemezed by stable thermal stratification with swell type 


Waves propagating through the area. 


rma Corporation System, Series 900 


eee stem Incorporation, Model 1054B 
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MONITOR 


SUPPRESSOR — AMPLIFIER 


Ww 






BRIDGE — 
AMPLIFIER 










(B) (from Davidson 1970) 


Figure 4. Temperature system (A) Diagram of the 
temperature probe, (B) Bridge-amplifier system 
Csonstant current) . 
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LV. STATISTICAL CONSIDERATIONS AND PROCEDURES 


A. SUMMARY OF JDF/CMF PROCEDURES 

Joint probability density function (JDF) and conditional 
Mean function (CMF) statistical procedures used in this study 
are based on descriptions by Holland (1968, 1972). The same 
procedures were applied by Frank (1971), Davidson and Frank 
(1973) and Thompson (1972) to small scale fluctuations over 
Waves. 

A JDF allows simultaneously observed pairs of values to 
be represented by contours on a two dimensional array. A 
trivariate statistical relationship is obtained by computing 
mmeomditLonal mean of a third variable as a function, CMF, 
of the two independent PHI) one vue ine fo gies) Di 
These statistics are determined for variables that are nor- 
memtized by their standard deviation. 

A JDF/CMF array appears in Figure 6A in which both the 
ordinate and abscissa are given in standard deviation (0) 
units. The array is centered on the mean value,. which is 
zero, of the two JDF variables. The bottom number on each 
eelienepresents the JDF (the probability x io of the joint 
feeunbpence or the pair). A nine point smoother was used to 
remove the perturbations viewed as being physically and 
Statistically insignificant. Statistical procedures neces~— 
Sary to obtain final JDF/CMF results are shown in aes d 
and will be described ‘peel ie in the following paragraph:. 


A complete discussion of the theory associated with the: 
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Figure 7. Block diagram of computer steps to 
produce a smooth JDF/CMF array. 


28 





procedures was presented by Holland (1968, 1972), and a more 
Bomplete discussion of its application in this study is pro- 


wieaged by Frank (1971). 


Bee BANDPASS FILTER 

The bandpass filter, shown as step 1, was an inverse 
transform filter. This was applied to isolate the signifi- 
cant frequency band (0.1 Hz to 0.8 HZ) associated with the 


ocean waves. Its response curve is shown in Figure 8. 





A n (Hz) ! 10 


(from Erank 1971) 


Figure 8. Response curve for inverse transform band- 
pass fiiter. 


C. NORMALIZATION OF DATA AND COMPUTATION OF STATISTICS 
Normalization of data, shown as step 2 in the flow dia- 


gram, was performed as described in subsection A. 
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D. JDF/CMF COMPUTATION 
The JDF/CMF array, shown as step 3 in Figure 7, is de- 
Semmuea@eiiebie Depanning of this Section. Procedures for 


this were described by Frank (1971) and Thompson (1972). 


Eee oLATISTICAL EDDY PERIOD COMPUTATION 

The "statistical eddy period" computed in step 3, was 
Suggested by Holland (1968, 1972) as an approach to reveal 
if the amplitude of the dependent variable such as tempera- 
ture varies nonlinearly with the amplitude of the independent 
variable such as wave height. Holland defined a "statistical 
eddy period" GA 27 a 
for any fluctuating meteorological variable y. oy is the 
Standard deviation of the meteorological variable and oy is 


meses agerivative. 


F. PHASE AMPLITUDE COMPUTATION 

Phase amplitude computation, shown as step 4 in the 
flow diagram, was determined from a JDF coordinate system 
with (n) as ordinate and its derivative (n) as the abscissa. 
This allowed presentation of the phase of the variable with 
respect to the wave. The phase angle was defined by an 
angle measured (by octants) counterclockwise (CCW) from the 
Norrzontal axeq in the (1,n) coordinate system. This is 
Shown in Figure 6B. A complete discussion of this may be 
Poumdmrm Holland (1968, 1972). 

The phase amplitude eee outlined in the above 


references were performed with respect to the waves fci che 
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Small Large 


Figure 9. Example of phase amplitude results for 
potential flow case. 





variables u, w, T and wT. Prior to summarizing by octants 
and plotting, the results were separated into wave amplitude 
masses defined by regiong on the JDF array. Three classes 
were considered, those being small amplitude (less than 1.5 
GQ), large amplitude (greater than 1.5¢) and the sum of the 
first two. This yielded, therefore, phase relationships 
with respect to small, large and all surface waves. 
Summarized phase amplitude results for each octant, as 
described in Figure 6B, are shown by horizontal bars that 
appear in Figure 9. Smooth lines were drawn through each 
set of horizontal bars such that areas above and below the 


resulting curves were equal. 


G. SUMMARY OF ANALYSIS PROCEDURES 

The previous discussion of JDF/CMF analysis procedures 
indicate that eine probability results are well suited for 
identification of fluctuations associated with waves. This 
type of result is particularly useful because it will reveal 
mom—Linear properties in fluctuations. These properties 
were alluded to by Volkov (1969) and Davidson and Frank 
no73). 

Poe ebemscmwown in the presentation of results, some of 
iemmost significant interpretations will be with respect to 
the phase Tit emidcom omit Oolcwam) Enese results do not repre- 
sent phase amplitude information as revealed by simple linear 
expansion of the fluctuations, i.e. spectral analysis, but 
iow mec ilemaectialenonlinear average fluctuations over differ- 


ent portions of the average wave. 


a 


2 





V. GENERAL WIND, WAVE, AND THERMAL CONDITIONS 


Data from four separate measurement periods were con- 
sidered. These periods were 19 May and 27 May and two on 
Eomiay L969, They will be referred to as periods 1 through 4 
respectively in the following discussion. The times for each 
period appear in Table 1 and Figure 5A along with information 
On the mean ten meter wind, the significant wave height, the 
meat flux (H) and Z/L. 

Wave histories in Figure 5A indicate that equilibrium 
conditions prevail with respect to the waves during all 
periods. Sea conditions were generally dominated by swell. 
Period 2 had the lowest significant wave height, 1.3 meters, 
While the other three periods ac Significant wave heights 
of about 2.6 meters. The wind speed at the 8 meter level 
ranged from 9.5 meters/second for period 4 to 11.8 meters/ 
second for period Dis 

The existence of equilibrium conditions with respect to 
atmospheric and sea conditions for all periods is evident, 
Pemitaps, by the fact that the relationship between oOT/T, and 
oe ueeaLe comparable to those observed by Volkov (1969). 
oly i. vs Sie results from those data are along the horizontal 
section of Yolkov's curye, Figure 10. Volkov (1969) has 
Suggested that this range of values corresponds to conditions 
with surge or nondeveloping sea. The curve shown in Figuré 


MWmismaereplot of Figure 2. 
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VI. PRESENTATION OF RESULTS 


Pec L ole eae So bATISOTECAL EDDY PERLOD" RESULTS 

Temperature spectra for periods 1, 3 and 4 appear in 
Paves log-log format of Figure ll. The spectral slope on 
these plots at high frequencies approaches a -5/3 value for 
all three periods. The -5/3 slope is expected for the 
inertial subrange of a passive scalar. 

Pairs of temperature, wl, and wave spectra for periods 
1, 3 and 4 are presented in linear-log profiles in Figure 12, 
Mmoeand 14. The left side of each Figure is the spectrum at 
the beginning of each period while the right side is the 
Spectrum occurring near the end. 

The most significant eHoeae reise tc of these spectra is 
the absence of any noticeable peak in the temperature spectra 
at the predominant frequency of the wave. Therefore, from 
an examination of spectra alone it might be mistakenly 
stated that ocean waves do not influence a passive scalar 
such as temperature. However, statistical results will be 
presented to show that temperature fluctuations are indeed 
influenced by wave motion below. 

Spectral results as presented above were not available 
Pormeperiod 2 . 

The "statistical eddy period" was computed as described 
in section 4 and the results are shown in Table 2. The 
Statistical eddy period for the temperature waveform is 


always less than that of the ocean wave. This strong]. 
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Suggests that potential flow theory will not accurately 


predict the air sea coupling in the boundary layer. 


B. JDF/CMF PRESENTATION FOR POTENTIAL FLOW CASE 

Ilempa Mec paiescomparLsons in the discussion of results 
will be with respect to potential flow theory. These 
Petrires are based on simple vertical displacement of a 
thermally stratified air column by the underlying progressive 
waves. If potential flow theory properly describes the 
coupling of the air column to the waves, then fluctuations 
i_meemperature (T), velocity (u), vertical motion (w) and 
heat flux (wT) observed by a fixed sensor would appear as 
shown in Table 1. JDF/CMF results would be those appearing 
in Figure 15 and phase amplitude results would be those 
appearing in Figure 9. 

A general description of the JDF patterns for potential 
flow case follows: 

1. Uniform circles; 90° phase shift between JDF variables 

meeipseseoOrientated 135°-315°, JDF variables are in 
phase 

Semeibipse, Oradentated 45°-215°, 180° phase shift in 
JDF variables. 

Dashed lines in all JDF/CMF results correspond to the 
Comoe sOrine Of the CMF. The orientation of various CMF 
Bee reais in each panel, as shown in Figure 15, is examined 


etoselyein the following interpretations. 
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3. 











eeeene Los 

Temperature and vertical motion are shown as the JDF 
and wave height as CMF in this panel. When temperature is 
maximum positive and wave height is maximum negative. Like- 
foc, the W maxima occur at every other inflection point of 
Ehe wave. 

27 ranel B 

In this panel the CMF zero isolines separates the 
positive and negative phase of the temperature waveforn. 

The JDF is the wave height and its derivative. The octants 
to be referred to were defined on page 30. This figure may 
be interpreted in a CCW direction starting in octant l. 
Temperature may then be interpreted to show that when the 
wave height is negative the temperature is positive and vice 
versa. “To state it another way, this panel represents a 
180° phase shift between the ocean wave and temperature wave 
forms. 

3. Panel C : 

The CMF zero isoline in this panel separates positive 
and negative vertical motion. The JDF is again determined 
from the wave height and its derivative. Interpretations of 
Sivompanel 1m 9a CCW direction from octant 1 indicates that 
if the wave height. is positive, w goes from positive to nega~ 
tive. This indicates a 90° phase shift between ocean waves 


and w wave forms during potential flow. 
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4, Panel D 
The u and T are presented in the JDF and the CMF is 
for wave height. In oectant 1 and 2 u and T are positive 
over the waves minima while in octant 5 and 6 u and T are 
negative over the wave maxima. 
5. Panel E 
Heat flux, wI, is the CMF and wave height and its 
derivative are the JDF in this panel. Beginning in octant 1, 
interpretation of this figure will describe heat flux in 
potential flow. During the positive phase of the wave wl 
completes a full negative and positive cycle with similar 
change for the negative phase of the wave. This indicates 
a 45° phase shift with respect to the wave. 
6. Panel F 
The CMF in this panel is u and the JDF is wave 
height and its derivative. Beginning in octant l, interpre- 
tation of the panel shows that when the phase of the wave is 
Mosttive the u waveform is negative and, therefore, a 180° 
phase shift. This is the same phase relationship shown for 
temperature in panel B. Attention will now be turned toward 
interpretation of results obtained from the four periods 


previously described. 


Cece RLOD I RESULTS 
Panels A and B of the JDF/CMF results, Figures 16, 1/7, 
Romande wltts be Considered in detail for the following 


Giscussion. A cursory examination of the JDF/CMF patter is 
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iepaneleA and Boot the above figures, in comparison to 

these in Figure 15 revealed that processes associated with 
potential flow are not the only ones that occurred. 
Specifically in Figure 16A there is a 135°-315° orientation 
Of the elliptical JDF pattern and the axis of the ellipse 
appears curvilinear. The elliptical pattern described above 
Occurred in all remaining subsets of this period, Figure l17A, 
18A and 19A, and during all periods considered in this study. 
This indicates approximately a 180° phase shift in w and T 
instead of the 90° observed in Figure 15 when potential flow 
was considered. The 90° CCW rotation of the CMF wave height 
pattern indicates that in general a w minima and T maxima 
Occur over the wave crest. 

The latter phase relationships are also apparent in the 
phase amplitude results in Figure 20. For potential flow 
temperature maxima occurred over the wave trough and w was 
90° out of phase with respect to the wave. A number of 
smaller perturbations appeared in the T waveform of Figure 
meee iney could be due to the turbulent contribution sug- 
gested by Volkov (1969). 

A possible physical interpretation of the phase rela:ion- 
ship noted above could be that at a fixed observation level 
the waves are causing a transfer of warmer air downward and 
colder air upward thereby inducing a net heat transfer down- 
Tanto ewas Met apparent, however, in the spectral 
results, Figure 12, 13 and 14, where it was observed that 


Oo eerecantwanectral peaks do not appear in or wT at 2 


46 





Mecdominate frequency of the wave spectrum. Temperature 
and w JDF and n CMF relationships throughout the remainder 
Geererdod | appear in panel A of Figures 17, 18 and 19. 
ihe CMF patterns in panel A of the above figures are no 
longer rotated 90° CCW from potential flow but are approxi- 
mately 45°. The 45° rotation is most evident in the large 
wave class region, greater than 1.50, of the JDF array. 
Figure 18A is an excellent example of this change in rotation. 
evaluation of the JDF/CMF results of Figure 18 in conjunc- 
tion with phase amplitude results in Figure 22 reveals that 
w and T are approximately 45° out of phase. Potential flow 
results had indicated a 90° phase shift. The temperature 
Maxima, considering only the larger crests, are almost in 
phase with the wave trough which indicates a tendency toward 
potential deena 

Changes in thermal stratification throughout the period 
may have been a factor in these results. The thermal strati- 
fication changed from near neutral to stable and returned to 
Meanmeneutral during the course of the period. The Z/L 
values changed from -0.03 to -0.20 and ended at -0.08. There- 
fore, with respect to stability and the above results one can 
Say that under greater stability there was a tendency toward 
potential flow as shown in Figure 18A and 22, and there were 
Fewer perturbations on the Beeuealay sinusoidal temperature 
waveform as shown in Figure 22. Furthermore, it was observed 


MmouiremroAwand 23 that as thermal stratification returned 
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toward neutral the JDF/CMF results deviate more from poten- 
tial flow conditions than was observed in Figure 18A and 22 
discussed above. 

Another feature apparent in panel A of Figures 16, 17 
and 18 is that the large wave class results near the outside 
of the array return to potential flow predictions more so 
than the small wave class results. 

Panel B of Figures 16, 17, 18 and 19 provided information 
on the phase relation of temperature with respect to the wave. 
Corresponding phase amplitude summaries for these JDF/CMF 
memays appear in Figures 20, 21, 22 and 23. JDF/CMF results 
in Panel B appear to be more sensitive to changes in thermal 
Stratification than those in panel A. For example, an examn- 
ination of Figure 16B reveals positive temperature over the 
wave crests while a similar comparison of Figure 1/7B reveals 
negative temperature over the wave troughs. Similarly an 
examination of Figure 18B reveals negative temperature over 
the wave crests as the period becomes more stable. ia 
general the JDF/CMF and phase amplitude patterns became less 
well defined as the period progressed and thermal stratifi- 
cation became more neutral. As an example, compare Figure 
isBawith 19B. 

Again the small wave class results differ more from 
potential flow prediction than the large wave class results. 

Weretideal action in eee tion with the ocean waves also 
appeacoulme panel G of Figures ieee eomand 19.  <~This series 


of panels cannot be readily interpreted for phase relet o2ns 
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between the w waveform and the ocean waves. Summarized phase 
amplitude results for w and the waves appear in Figures 20, 
21, 22 and 23. The results of these figures also show that 
there is a W maxima anaes of the wave crest. For the poten- 
tial flow case, there would be a 90° lead. 

Terereciuleswbor the Joint occurrence of u and T with 
respect to the CMF of the wave height appear in panel D of 
Figures 16, 17, 18 and 19. The 045°-215° orientation of the 
JDF ellipse indicates u and T are positively correlated as 
predicted by potential flow. However, the curvilinear shape 
of the axis of the ellipse represents a deviation from 
Peienitial Plowwemeory. Lhe CMF zero isoline for wave height 
ieepanel 1 is rotated 45° CW from the potential flow orienta- 
tion. This indicates a phase shift of both u and T with 
respect to the wave. Summarized phase amplitude results for 


Meand <I appear in Figures 20, 21, 22 and 23. 


eee Ee RIOD 2 RESULTS 

Thermal stratification was slightly stable at the begin- 
ning of the period and became near neutral as the period 
progressed as indicated by Z/L values changing from -0.0- to 
Be eee Ooulystne tirst two subsets of the form will be dis- 
cussed for this period because all subsets, in contrast to 
BiGacm juste described for period 1, are quite similar. Sub- 
set 3 and 4 results appear in Appendix A in Figures 32, 33, 
34 an Bs 

In general results reveal deviations from potential low 


conditions. This can be readily noted by comparing pati rns 
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(A) (B) (C) 
Figure 21. Phase amplitude results for period 1 subset 2 
(a) small waves, (b) large waves, and (c) all waves. 
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Pe tire ower nase amplitude results for period 1 subset el 
fapmectiviiewavesc.) (b) large waves, and (c) all waves. 
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Fagure 23. Phase amplitude results for period l subset 4 
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Per ieures 24 and 25 with results in Figure 15. eee 
eaeivy, the CME zero isoline in Figure 24A has a 45° CCW rota- 
monet rom potential flow orientation, Figure 15A. A similar 
rotation was observed in Figure 18A which corresponds to the 
more stable conditions in Period ll. 

The T and w/JDF results are elliptical in Figure 24A as 
in Period 1 and has a 315°-135° orientation. This again 
Suggests a 180° phase shift between T and w in contrast to 
ae 90° predicted by potential flow theory shown in Figure 
15. 

Relationships between features in this period and the 
thermal stratification can be further examined by a broad 
comparison of similar JDF/CMF panels in Figure 24 and 25. 
The Pris in Figure 24 appear to have fewer small scale 
perturbations than those in Figure 25. Phase results for 
this period eee tied in the phase amplitude results in 
Figures 26 and 27. It can be noted that the results in Fig- 
Meee27 in contrast to the JDF/CMF results in Figures 25 are 
more readily examined even though small scale perturbations 
occur. An example is the temperature waveform in Figure 2/7. 
inmermal stratification was near neutral during the second 
subset corresponding to Figure 27. This deviation from a 
sinusoidal waveform, as for Period 1, is considered to be 
ine to turbulence. 

With respect to different wave classes, small wave 
classes appear to approach potential flow features more than 


Pte latrermmwave classes. This feature is shown by the’ ve 
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(emote itezeroOmicoline in Figure 24A. This is in contrast 
to their tendency to deviate from potential flow as was ob- 
served in Period 1 during more thermally stable conditions. 
Mencacdremenatedid not appear in Period 1 but appears in 
this period is the tendency for the temperature waveform to 


have sharper crests and troughs than the ocean waves. 


Ewe FERIOD 3 RESULTS 

Hydrostatic conditions throughout this period were stable 
Peed 2/4 value of approximately -0.12. Results for this 
period were quite similar to those presented for subsets 2 
and 3 of period 1, Figures 17 and 18. The latter subsets 
also occurred during more stable conditions. 

Peccmatmcumec subsets of this perlod were similar, 
only subset 1 will be discussed in detail. Results for sub- 
sets 2 and 3 are in Appendix A as Figures 36,37,38 and 39. 

The JDF of u and T and the CMF of wave height appear in 
Paoure 260A. <The CMF zero isoline in Figure 28A is rotated 
approximately 70° CCW from potwentialeflow results, Figure 15. 
Phase amplitude results in Figure 29 indicate the above 
rotation is due to the temperature waveform being shifted 
approximately 45° frometne wave s trough. This is to be 
expected in view of previous results when thermal stability, 
ieeerms Of Z2Z/l, values, was considered. During period l, 
wy levalues varied from -0.03 to -0.20 to -0.08 and maxima in 
the temperature phase amplitude varied from 180° out of phase 


with the wave trough, Figure 23, to approximately in ph se 
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Figure 25. JDF/CMF Results for period 2 subset 2. 


61 








4. 
i) 
Sa) 





/sec 


mcm 
Sa) 











ae 
> 
© 


/sec 





ts Ci 








aSke 
a a 
wv Sex 
w 
a A 
j } 
\ 7a 
eee rie 
(A) (B) (C) 
Figure 26. Phase amplitude results for period 2 subset l 


(a) small waves, (b) large waves, and (c) all waves. 
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Figure 27. Phase amplitude results for period 2 
Subsets (a) small waves, (b) large waves, and (c) 
all waves. 
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Peuiethe wave LEough, Figure 20. Period 3,however, has a 
Pomscictentu2;/ lb evalwe of about -0.12. Therefore, phase dif- 
ferences between temperature and waves during period 3 were 
between the extremes noted in period 1. This can be readily 
observed by a comparison of the phase amplitude results for 
T and waves in Figure 20 and 23 to those in Figure 29. 

With respect to wave amplitude classes, small and large 
waves appear to have the same deviation from potential flow. 
This is most evident by a comparison of the wave height CMF 
Metro isoline in Figure 28A and in Figure 15. Recall that 
during the changing stability associated with period l, only 
the small wave amplitude classes had a tendency toward 
metential flow, Figure 18A. This is in contrast to the 
results during the near neutral conditions during period 2 
where small waves had a tendency toward potential flow, 
Figure 24A. 

Panel 28B and 28C results show the CMF variables T and w 
respectively compared to the JDF of the wave height and its 
derivative. CMF zero isolines in both panels are shifted 
ate wrronepotential flow, Figure 15. In addition, panels 
Sooeand) © display fewer small perturbation than similar re- 
Pieomobtaimed in period 1, Panel B or C of Figures 16, l/, 
18 and 19 or period 2, Panel B and C of Figures 24 and 25. 
tis sacain was probably associated with the turbulence occur- 


mene during the period. 
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Figure 28. JDF/CMF Results for period 3 subset 1 (a) 
small waves, (b) large waves, and (c) all waves. 
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Peure sz. hase amplitude results for period 3 subset 1 
(a) small waves, (b) large waves, and (c) all waves. 
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Pe PERIOD 4 RESULTS. 

Results are similar to those in Period 3. The value of 
Mie thermal stabrlity parameter, Z/L, averaged -0.12 in this 
period but approaches neutral conditions near the end of the 
period. 

The period ended with a Z/L value of -0.05. Even though 
the stability decreased during period 4 the results were 
quite similar to those obtained for period 3. However, the 
near neutral conditions did not occur until the very end of 
the period and, thereby, suggest that some time lag occurred 
after the long interval of stability in period 3 before the 
wind wave coupling adjusted to the new environment and the 
analytical method used could detect charges. Due to the 
Semitarity of results for this period, Only ubne first of. the 
3 subsets will be used for discussion. Results for the last 


Buon subsets are presented in Appendix A, Figure 40, 41, 42 


and 43. 

Pecuilts fon this period, shown in Figure 30, are similar 
MomEiiose for period 3. For example, in Figure 30A, the 
temperature maxima is 45° from the wave's trough as was the 
case during period 3, Figure 28A. Phase amplitude results 


wee shown in Figure 31. 
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Poqimemov0ueJbr/CME Results for period 4 subset l. 
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Pupiner ote Phase amplitude results for period 4 subset l 
(a) small waves, (b) large waves, and (c) all waves. 
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Vil es OONCLUSIONS 


iicmcoeucyeids Deen an investigation of turbulence data 
obtained during near neutral to stable conditions. The tenm- 
perature fluctuations, were described with respect to the 
wave's influence over natural ocean waves. The statistical 
procedure, JDF/CMF, was used which made possible interpreta- 
tion of nonlinear properties of wave induced temperature 
mluctuation. 

The general characteristic of observed temperature wave- 

forms leads to the following conclusions. 

1. The most important property shown is that waves 
influence a scalar field such as temperature as high 
as 8 meters. Therefore, observational studies con- 
ducted eee ocean waves must consider the wave's 
influence on a scalar field. 

2. The suggestion by Volkov (1969) that a temperature 
observation at a fixed level above a wave surface is 
composed of a mean, a fairly well organized fluctua- 
buonmmawentOo waves, and a fluctuation due to turbu- 


lence appears to be valid. All periods show small 


fluctuations superimposed on a more uniform sinusoidal 


wave induced fluctuation. 

3. The phase shift of temperature away from wave crests 
is 0-135° with an average about 84° for all periods. 
This is in contrast to the 180° phase shift pimemicted 


bye potential flow. 
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Tee tomtctiperatwre profile leads the vertical velocity 
although it may have a maximum of 45° on either side 
of the wave crest. This is in contrast to the 90° 
shift predicted by potential flow. 

5. As stated in the background, Stewart (1969) did not 
consider waves capable of influencing a scalar tem- 
BewamuGemitroid. ~ therefore, he considered vertical 
velocity and T phase angles near zero or 180° indica- 
tive of mostly turbulence and angles + 90° indicative 
of waves. This study has shown that under stable 
conditions neither extreme is attained. The reason 
for not consistently being at either extreme is 
possibly due to the ee Piece defile twat 2Om,. 

6. It also appears for sie one Study e@tmat if the 
atmosphere is neutral or very stable for an extended 
period, sea~-air coupling will tend toward, but 
probably Ree ek. moLentigie t low characteristic. 
In between these two extremes, the greatest deviation 
Eprom potential flow will be seen. 

iiicmsottdyeshows that JDF/CMF statistical technique is 

Mmeoeyepromising and should be applied to other data obtained 

Strict mecnsame period. Such a study would determine if the 
conclusions of this study are consistent. Unstable periods 

should be included if available. 

One possible application presented for this type of 

Baldev aiomatomice in predicting laser propagation. With 


DPeeemortdidcscmonmenis type it should become possible t¢ 


ay 





MeciiereENcmmciperatinre fluctuations under different thermal 
emir bcateomse) taererore, with a knowledge of the thermal 
stratification and the wave field below, a more accurate 

refractive index corrected for the wave's influence could be 


provided to laser users. 
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JDF/CMF Results for period 2 subset 3. 
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Figure 34. Phase amplitude results for perrod 2 
subset 3 (a) small waves, (b) large waves, and 
(c) all waves. 
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Figure 35. Phase amplitude results for period 2 
subset 4 (a) small waves, (b) large waves, and 
(c) all waves. 
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Wien becults for period 3 subset 2. 
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Peete Dr /CMP Results for period 3 subset 3. 
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Exegure 38. Phase Miplimimdemrectlts £or period 3 subset 2 
(a) small waves, (b) large waves, and (c) all waves. 
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Figure 39. Phase amplitude results for period 3 subset 3 
(a) small waves, (b) large waves, and (c) all waves. 
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Figure 40. 
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Peter tO /CME Results for period 4 subset 3. 
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Premre 42. Phase amplitude results for period 4 subset 2 


fmecmallewaves. (€b) large waves, and (c) all waves. 
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(a) small waves, (b) large waves, and (c) all waves. 
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